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Abstract. We study a one-dimensional reaction-diffusion system which describes an 
isothermal autocatalytic chemical reaction involving both a quadratic {A + B —>■ 2B) 
and a cubic {A + 2B — > 3-B) autocatalysis. The parameters of this system are the ratio 
D = Db/Da of the diffusion constants of the reactant A and the autocatalyst S, and 
the relative activity k of the cubic reaction. First, for all values of D > and /c > 0, we 
prove the existence of a family of propagating fronts (or travelling waves) describing the 
advance of the reaction. In particular, in the quadratic case /c = 0, we recover the results 
of Billingham and Needham [BN]. Then, if D is close to 1 and k is sufficiently small, 
we prove using energy functionals that these propagating fronts are stable against small 
perturbations in exponentially weighted Sobolev spaces. This extends to our system 
part of the stability resiilts which are known for the scalar Fisher equation. 
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1. Introduction 



We consider the reaction-diffusion system 

dtu{x,t) = d'^u{x,t) — u{x,t)v{x,t) — ku{x,t)v{x,tY , 
dtv{x,t) = Ddlv{x,t) + u{x,t)v{x,t) + ku{x,t)v{x,t)'^ , 



(1.1) 



where u,v are nonnegative functions of (x^t) G IR x IR_|_, and D > 0, k > are 
constant parameters. This system describes (in dimcnsionless variables) an isothermal 
autocatalytic chemical reaction of mixed order, involving both a quadratic {A+B 2B) 
and a cubic {A + 2B — 3i?) auto catalysis, see [HPSS], [BN]. Here, u and v are the 
concentrations of the reactant A and the autocatalyst S, and D = Db/ Da is the ratio 
of the diffusion constants. The parameter k measures the contribution of the cubic 
autocatalysis to the whole reaction. Of particular interest are the purely quadratic case 
/c = 0, and the purely cubic case "/c +cxd" which corresponds, after a rescaling, to 
Eq.(l.l) with k = 1 and without quadratic terms. 

The dynamics of the system (1.1) on a bounded domain f2 C H with homogeneous 
boundary conditions is well understood [Ma], [HY]. For any initial data uq^vq G L°°(f2), 
the solution {u{t)^v{t)) = {u{-,t),v{-,t)) G L°°(f])^ is uniformly bounded for all times 
and converges as t — > +oo to a uniform steady state {u*,v*) satisfying u*v* = 0. On the 
other hand, if uo,vo G L-^(]R) nL°°(]R), the solution {u{t),v{t)) of the system (1.1) on 
the whole real line stays bounded for all times [BKX] and converges uniformly to zero as 
t — > +00. In the purely cubic case, a very detailed description of this convergence can be 
found in [BX], [BKX]. Finally, if uo,vo G L°°(IR) only, then the solution {u{t),v{t)) exists 
for all times, and stays uniformly bounded if D < 1 [MP]; if D > 1, uniform boundedness 
is an open problem, but an upper bound is known which diverges extremely slowly as 
t — > +00 [CX]. Of course, very little is known about the behavior of the solutions in 
this general situation. 

In this paper, we investigate the existence and stability of propagating fronts (or 
travelling waves) for the system (1.1). These are uniformly translating solutions connect- 
ing the stable steady state («, v) = (0, 1) at a; = — oo to the unstable state (m, v) = (1, 0) 
at a; = +oo. Thus we look for solutions of (1.1) of the form u{x,t) = a{x — ct), 
v{x,t) — I3{x — ct), where c > is the velocity of the front. The nonnegative functions 
a, (3 satisfy the system 



a" {x) + ca' {x) — a[3 ~ ka0^ = , 
D(3"{x) + c/3'{x) + a/3 + ka0^ = , 



(1.2) 
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together with the boundary conditions 

(a(-oo),/?(-oo)) = (0,1) , (a(+oo),/9(+oo)) = (1,0) . (1.3) 

Remark. It is easy to verify that the nonnegative time-independent solutions of (1.1) 
are exactly the uniform steady states (-u, v) = (a, b) with a > 0, 6 > satisfying ab = 0. 
Moreover, a necessary condition for the existence of a trajectory of (1.2) connecting 
(a_, b-) at X = — oo to (a+, b+) at x = +oo is that a_ + 6_ = a+ + b+. Indeed, adding 
the two equations in (1.2) and integrating with respect to x, we obtain the conservation 
law 

a'{x) + ca{x) + D(3'{x) + c(3{x) = const. , (1-4) 

and the assertion follows by taking the limits x — > ±oo. Therefore, any heteroclinic 
orbit of (1.2) must connect (0, a) to (a, 0) for some a > 0. Now, the system (1.1) is 
invariant under the scaling transformation 

u{x,t) ^ X^uiXx.X^t) , v{x,t) ^ X^v{Xx,X^t) , k^k/\^, (1.5) 

for all A > 0, so there is no loss of generality in assuming that a = 1. This explains the 
choice of the boundary conditions (1.3). 

Existence of solutions to (1.2), (1.3) has been studied by Billingham and Needham 
[BN] for all D > in the cases A; = and k = +oo. First, they show that any solution 
satisfies the following properties for all x eJR: 

0<a{x)<l, a'{x)>0, j a{x) + (3{x) < 1 if < 1, 

0</3(x)<l, P'{x)<0, \a{x) + P{x)>l ifD>l. 

Then, in the purely quadratic case, they prove that a travelling wave exists if and only 
if c > 2\/D, and is unique up to a translation in the variable x. Finally, in the purely 
cubic case, they argue that a propagating front exists if and only if c > V2{D), where 
V2{D) is some increasing function of D satisfying f|(l) = l/\/2, V2{D) — 0{D) as 
D ^ 0, and f|(-D) = 0{VD) as D ^ +oo. In this latter case, their argument relies in 
part on numerical calculations. 

On the other hand, if D = 1, the existence of travelling waves has been proved for 
all A; > 0. Indeed, it follows from (1.6) that a{x) + f3{x) = 1 in this case, so that (3 
satisfies the single equation 



13" {x) + cl3\x) + I3{x)(l - I3{x)){l + kl3{x)) = , 



(1.7) 
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together with the boundary conditions oo) = 1, /3(+oo) = 0. This problem can be 
studied by usual phase space techniques, and is known to have a nonnegative solution 
if and only if c > c*(/c), where 



see [BBDKL], [vS]. If A; < 2, the front with minimal speed is often called "pulled" or 
"linear", because its velocity c* and its decay rate at infinity can be determined from 
the linearized equation ahead of the front. In the converse case, it is called "pushed" or 
"nonlinear" . 

In the general case D > 0, k > 0, the situation seems more complicated, and our 
results are still incomplete. If D > 1, we can still prove the existence of a minimal 
propagation speed for the travelling waves: 

Theorem 1.1. Let D > 1, k > 0. Then there exists c* = c*{D,k) > such that 
Eqs.(1.2), (1.3) have a nonnegative solution if and only if c > c*. This solution is 
unique up to a translation and satisfies (1-6). Moreover, c*{D,k) is a non-decreasing 
function of k and satisfies 



where c{D, k) is given by Eq.(2.13) below. In particular, c{D,k) < VD{Vk + ^/TJk) 
for allD>l and all k > (3D - 1)/(3D - 2). 

An immediate consequence of Theorem 1.1 is: 

Corollary 1.2. Let D > 1. There exists k* = k*{D) such that the minimal speed 
c*{D,k) defined in Theorem 1.1 satisfies c* = 2VD if k < k*(D) and c* > 2VD if 
k > k*(D). Moreover, one has 



for all D>1. 

This result says that the front with minimal speed is "linear" if k < k* and "non- 
linear" if k > k*. A numerical determination of the curve k*(D) is shown in Fig. 1 
below. 




(1.8) 




(1.9) 
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On the other hand, if L> < 1, we do not know whether the set of values of c for 
which a propagating front exists is always an interval of the form [c*, oo). However, for 
all D < 1 and all A; > 0, we do know that a nonnegative solution to Eqs.(1.2), (1.3) 
exists if c is sufficiently large and does not exist if c < 2\/D: 

Theorem 1.3. Let Q < D < 1, k > Q. Then Eqs.(1.2), (1.3) have a nonnegative 
solution if 

^ f 2y/D ifk < 2 , 

^ ~ \ ^/D{y/k/2+ y/2/k) ifk>2. 

This solution is unique up to a translation and satisfies (1.6). Conversely, there exists 
no nonnegative solution of Eqs.(1.2), (1.3) if c < 2\/D. 

Remark. In particular. Theorem 1.1 and Theorem 1.3 show that, D > 1 and k < 
(3-D — 1)/ {3D — 2), or if D < 1 and k < 2, a nonnegative propagating front exists if and 
only if c > 2^JD. This extends the result obtained by Billingham and Needham [BN] 
for /c = 0. 
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Fig. 1: The curve k*{D) separating the parameter regions where c*{D, k) = (hnear fronts) and 

c*{D, k) > 2-v/D (nonUnear fronts). The existence of this curve is asserted by Corollary 1.2 for D > 1 
only, but numerically this separatrix can be observed for all D > 0. Note that k*{D) > 2 if < D < 1, 
in agreement with Theorem 1.3. 



It follows from Theorem 1.1 and Theorem 1.3 that, for all -D > and all /c > 0, 
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the system (1.1) has a one-parameter family of uniformly translating front solutions, 
indexed by the velocity c. A natural question to address is whether these propagating 
fronts are stable against sufficiently small perturbations in appropriate function spaces. 
Again, the case D = 1 is easier and can be treated separately. Indeed, if D = 1, the 
system (1.1) can be written as 

d+wix, t) = dlwix, t) , 

(1.10) 

dtv{x,t) — dlv{x,t) + v{x,t){w{x,t) — v{x,t)){l + kv{x,t)) , 

where w{x,t) = u{x,t) + v{x,t). In these new variables, the propagating fronts are 
given by w{x,t) = 1, v{x,t) = P{x — ct), where /? is the solution of (1.7). Therefore, if 
we consider initial data of the form wo = l + f,vo = P + g, where /, g are sufficiently 
localized perturbations, then the solution of (1.10) will satisfy w(x, t) = l + 0(t~^^'^) as 
t — > +00, so that the behavior of v{x, t) for large times will be governed by the nonlinear 
diffusion equation 

dtv{x,t) = dlv{x,t) + v{x,t){l - v{x,t)){l + kv{x,t)) , (1.11) 

up to a remainder which can be controlled rigorously [Fol]. Now, the stability of the 
travelling wave solutions of Eq.(l.ll) has been intensively studied by many authors 
[AW], [Sa], [Ki], [EW], [BK], [Ga2]. In particular, for all /c > and all c > c*(/c), each 
individual front is known to be asymptotically stable against perturbations which decay 
to zero sufficiently fast — at least as fast as the front itself — as x ^ +00. The decay 
rate in time of the perturbations is polynomial or exponential depending on the choice 
of the function space. In addition, if A; > 2 and c = c* (the "pushed" case), the family 
of all translates of the front is orbit ally stable against perturbations which decay even 
slower than the front itself as x — > +00. 

In the general case -D 7^ 1, the reduction to a single equation is no longer possible, 
and much less is known about the stability of propagating fronts. We shall restrict 
ourselves in the sequel to the situation where D is close to 1 and k is close to 0, see the 
remarks after Theorem 1.4 below for a discussion of these limitations. In this parameter 
region, we know from Theorem 1.1, Theorem 1.3 that a propagating front (a,/?) exists 
if and only if c > 2\/D. Since a{x) — > and /3{x) — > 1 as x ^ —00, we shall assume 
(without loss of generality) that P{x) — a{x) > 3/4 for all a; < 0. Setting 



u{x, t) = a{x — ct) + f{x — ct, t) , v{x, t) = P{x — ct) + g{x — ct, t) , 
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and inserting into (1.1), we obtain the evolution equations for the perturbation {f,g) 
in the moving frame 

dtf = dlf + cd^f - /3(1 + A;/3)/ - a(l + 2kp)g - N{f, g) , 
dtg = Ddlg + ca,^ + a{l + 2k(3)g + (3{l + /c/3)/ + iV(/, ^) , 

where iV(/, g) = fg + k{2l3fg + ag^ + fg2y 

As in the scalar case, it is necessary to use weighted spaces which force the per- 
turbations (/, (7) to decay to zero sufficiently fast as a; ^ +00 [Sa]. For any s > 0, we 
consider the Hilbert spaces Xg, Yg of real functions on IR defined by the norms 

f \h{x)\\l + e'ndx, \\h\\l^ = ll/^lli + ll/^'lli , (1-13) 

^IR 

where ' denotes the (space) derivative. We also note X'^ = Xs x Xg, =^5X1^3. 
Then a direct calculation shows that the origin (/,(/) = (0, 0) in (1.12) is linearly stable 
in X^ only if Ds^ — cs + 1 < 0, see [Gal], [GR] for a similar discussion. Therefore, we 
must choose s = l/\fD if c = 2\fD. If c > 2\/D, s can be chosen in a whole interval, 
but the biggest perturbation space corresponds to the choice 

s = ^{c-^c^-AD), (1.14) 

which we shall always assume in the sequel. Note that this value corresponds to the 
exponential decay rate of both a{x), (3{x) as x ^ +00. 

With these definitions, we can state our main result: 

Theorem 1.4. There exist do > and ko > such that, for all D G [1 — (io, 1 + 
do], aii k G [0, /cq] and all c > 2^/D, there exist cq > and Kq > 1 such that the 
following holds: for all (fo^go) G Y^ satisfying \\{foj go)\\Y^ < eo, where s is given by 
(1.14), Eq.(1.12) has a unique global solution (/, g) G C°([0, +00), Y^)nC\{0, +00), Y^) 
satisfying (/(0),gf(0)) = {fo,go)- Moreover, one has 

\\{fit),gm\Yi < i^o||(/o,5o)|ln^ , (1.15) 

for all t >0, and 

hm \\{dJit),dMt))\\x2 = 0. (1.16) 

t— »-|-oo " 
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Remarks. 

a) In particular, Eqs.(1.15), (1.16) imply that the perturbation {f{t),g{t)) of the front 
converges to zero uniformly in the following sense 



but Theorem 1.4 does not give any pointwise decay rate in time for the perturbations. 
However, our proof could be extended to provide some (non-optimal) decay rate at the 
expense of using higher order energy functionals, see for example [FS]. On the other 
hand, a detailed study of the linearized equation (1.12) (with N{f,g) = 0) including a 
careful determination of the optimal decay rate in this case can be found in [Fo2] . 

c) A important open problem is whether the result of Theorem 1.4 can be extended to 
other values of the parameters -D, /c, using possibly different perturbation spaces. The 
difficulty we encounter when D is far from 1 is related to the "Turing phenomenon" in 
the theory of pattern formation [Tu], which shows how a stable equilibrium point of a 
reaction system can be destabilized by diffusion if the components in the system have 
very different diffusion rates. The question is therefore whether this mechanism actually 
leads to instabilities in our system if either D <^ 1 or D » 1. On the other hand, the 
difficulty we have when A; » 1 is more technical in nature: the spectral analysis of 
the linearized operator in (1.12) becomes difficult when k is large, because we have 
to preclude the existence of unstable eigenvalues. This problem already exists when 
D = 1, but in this case the reduction to a single equation allows one to use some results 
from the theory of Schrodinger operators which do not extend to systems. Another 
interesting question is therefore whether Theorem 1.4 holds true for all A; > 0, at least 
if D is close to 1, and whether the orbital stability result for the family of translates of 
the "pushed front", which is known ior D = 1, k > 2, c = c*{k), has any analogue in 
the general case D ^ 1. 

The rest of this paper is organized as follows. In section 2, we prove the existence 
of propagating fronts (Theorem 1.1 and Theorem 1.3) using phase space techniques. 
In section 3, we construct energy functionals which allow us to show that these front 
solutions are stable against sufficiently small perturbations in (Theorem 1.4). 



hm sup{\f{x,t)\ + \g{x,t)\){l + e'^) = 0. 



b) In addition to (1.16), the proof will show that 
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2. Existence of Propagating Fronts 

In this section, we study the existence of nonnegative solutions to the system (1.2) 
satisfying the boundary conditions (1.3). Our method, which follows closely Billingham 
and Needham [BN] , relies on the construction of invariant regions in a three-dimensional 
phase space. When D < 1, A; > 0, we show that a propagating front exists if c is large 
enough, and does not exist if c is too small, thus proving Theorem 1.3. A similar result 
holds when D > 1, but in addition we can show that the existence of a propagating front 
for c > 0, A; > implies the existence for all c' > c. A;' < k, thus proving Theorem 1.1. 

We first note that any solution of (1.2), (1.3) corresponds to a heteroclinic orbit of 
a three-dimensional dynamical system, see [BN]. Indeed, setting (3' = w and using the 
conservation law (1.4) to eliminate a\ we see that a,(3,w satisfy the equations 

a' = c{l-a- p) - Dw , P' = w , w' = -^{cw + aP(l + k^)) , (2.1) 

together with the boundary conditions {a, /3, w){—oo) = Pi = (0, 1, 0), {a, (3, w){-\-oo) = 
P2 = (1,0,0). 

Straightforward calculations show that the fixed point Pi is a saddle for the dy- 
namical system (2.1), with one positive and two negative eigenvalues given by 



C2 



X± = --±^- + l + k, Ao = --. (2.2) 



The eigenvector corresponding to A_|_ > is 



On the other hand, the fixed point P2 is a sink, with three negative eigenvalues given 

by 

= ^(-c ± \/c2 - 4L>) , iio = -c. (2.4) 

Note that iJ,± > fj,o ii D > 1/2 and 1 + {D — l)c^ > 0. The eigenvectors corresponding 
to ii±, jJLQ are respectively 

w± = f-^±^E± , 1, , wo^ (1,0,0) . (2.5) 



Let V be the one-dimensional global unstable manifold of Pi , and let V+ C V be the 
invariant manifold which coincides with V n Q+ in a small neighborhood of Pi, where 
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= {{a, P^w) \ a > 0, f3 > 0}. In view of the preceding remarks, proving the existence 
of a nonnegative propagating front amounts to showing that V+ C Q+ fl Q{P2), where 
$7(^2) is the basin of attraction of P2- Following [BN], we shall prove this, for some 
values of the parameters D, c, k, by constructing an invariant region TZ C Q+ H 0(^2) 
such that V+ C TZ. Since everything is known when D = 1, we shall only consider the 
cases < D < 1 and D > 1. 



2.1. The Case D <1 

Let D e (0, 1). Given A> 0, we define the (closed) region TZi C Q+ by 

7^l = {(«,/?,«;) |0 < (3<1 , 0<a<l-P , -AP{1- (3) < «; < 0} . (2.6) 
We also note 7^l = \ {Pi}. 

Lemma 2.1. Assume that the parameters D<l,c>0, k>0, A>0 satisfy 

DA^ - cyl + 1 < , k- 2DA^ < . (2.7) 
Then TZi is invariant under the Sow of (2.1 ), and V+ C TZi C ^(^2)- 

Proof. To prove the invariance of TZi, we show that the vector field (2.1) at any point 
P G dTZi is directed into TZi or is parallel to the surface of TZi. This is easy to verify 
for the faces a = 0, a = 1 — P, and w — 0. On the last face w = —AP{1 — P), we have 

w' + AP'{1 -2P) ^ w(^+ A{1 - 2/3)^ - ^aP{l + kp) 

= APil -P){^- Ail - 2/3)) - la/3(l + kp) (2-8) 

> ^P{1 - p) {-{DA^ - + 1) + p(2DA^ - k)) , 

since a < 1 — p. Obviously, the right-hand side of (2.8) is nonnegative if (2.7) holds, 
hence TZt is invariant under the flow of (2.1). 

Due to this invariance, to prove that V+ C TZi it suffices to verify the inclusion in a 
small neighborhood A/" of Pi. Since V+ fl A/" C {Pi + ev+ + 0{e^) | < e < eo} for some 
eo > 0, it follows from (2.3), (2.6) that V+ nTV C Tei if 

^±£h± < 1 , and X+ < A. (2.9) 
c + X+ 
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The first inequahty is always satisfied since D < 1, A-(- > 0, and the second one follows 
from (2.2), (2.7). Indeed, adding the two inequalities (2.7), we obtain DA^ + cA — (1 + 
k) > 0, while A^. + cX+ - {1 + k) =0 by (2.2). Since D < 1, this implies A > A+. 
Therefore, we have shown that V+ C TZi, hence V+ C TZi since Pi ^ V+. 

Finally, let P e T^-i, and let 7(0;) = {a, P,w){x) be the solution of (2.1) satisfying 
7(0) = P. Since 7(0;) e TZi for all a; > 0, we have a'{x) > 0, p'{x) < for all x > 0, 
hence a{x) a e[0,l], f3{x) ^ ^5 e [0, 1) as a; ^ +00 (note that ^ < /?(0) < 1.) From 
(2.1), we deduce that w{x) w = —c~^aP{l + kp) as a; — > +00. Now, the fixed point 
(a, w) has to satisfy w = 0, aP = and a + /3 = 1, hence a = 1, P = 0. This proves 
that ^{x) — > P2 as a; — > +00, hence TZi C 0,{P2). □ 

Proof of Theorem 1.3. From Lemma 2.1, we know that a nonnegative propagating 
front exists if the conditions (2.7) can be satisfied for some ^ > 0. This front is unique 
up to a translation in x since the unstable manifold V+ is one-dimensional, and the 
properties (1.6) follow from (2.6) or can be verified directly as in [BN]. Now, if A; < 2, we 
choose A = 1/VD, and (2.7) holds for all c > 2y/D. li k > 2, we choose A = y/k/{2D), 
and (2.7) holds for all c > ^/D{y/k/2+ y/2/k). This proves the first part of the result. 

Conversely, if c < 2\/D, the eigenvalues /U± in (2.4) become complex, and it is easy 
to show that no trajectory of (2.1) can stay in and converge to P2 as a; — +00, except 
on the invariant line P = w = which does not intersect V+. Therefore, no nonnegative 
front solution can exist if c < 2\fD. This concludes the proof of Theorem 1.3. □ 

Remark. Following the same lines, one verifies that the region 

Tlx = {{a,P,w)\^<P< 1 , {1-P){\-EP) <«< 1-/3 , -AP{\-P) < ti; < 0} 

satisfies the conclusion of Lemma 2.1 if £^ > and if the additional condition ^(1 — 
D) — (c — A^)E < is fulfilled (see also Lemma 2.2 below.) In particular, if D > 1/2, 
c > 2^fD and /c < 2, we can choose A = 1/VD and E = {1 - D)/{2D - 1). This shows 
that the propagating front satisfies the bound 

(1 - /5) (1 - <a<l-p. (2.10) 
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2.2. The Case D > 1 : Bounds on the Critical Speed 

Let i:) > 1, ^ > 0, E > 0. We define the region 7^2 C Q+ by 

7^2 = {(a, tu) I < /? < 1 , 1-P<a<(l-P)(l + Ep) , -AP{1 -P)<w<0} . 
We also note 7^2 = Tea \ {Pi}. 

Lemma 2.2. Assume that the parameters D>l,c>0, k>0, A>0, E>0 satisfy 
DA^-cA + l<0, E{l + k) + k-2DA^ < , A{D -1) - E{c- A) < . {2A1) 
Then TZ2 is invariant under the Row of (2.1), and V+ C 7^2 C fl{P2)- 

Proof. We proceed as in the proof of Lemma 2.1. First, it is easy to verify that 
the vector field (2.1) is directed into TZ2 on the faces a = 1 — P and w = 0. When 
w = —AP{1 — P), we have as in (2.8) 

w' + Ap'il-2f3) = AP(1 -P){-^- ^(1 - 2/?)) - ^ap(l + k(3) 

> ^(3{1 - P) {cA - DA\1 - 2/3) - (1 + EP){1 + kp)) 

> ^P{1 - P) {-{DA^ -cA+1) + p{2DA^ -k-E{l + k))) > , 

since a < {1 - p){l + Ep) and P'^ < p. Finally, when a={l- p){l + Ep), we find 

a' + P'{1 -E + 2EP) = -cEP{l -P)+w{l-D-E + 2EP) 

< P{1 - P){-cE + A{D-1 + E))<{), 

since —Ap{l — P) < w < 0. This proves that 7^.2 is invariant under the flow of (2.1). 
To show that V+ C 7^2, it is sufficient to verify that 

1 < ^ ^ < 1 + E , and X+ < A , (2.12), 
c+ A+ 

see (2.9). The first inequality is obvious since D > 1. To prove that A-|- < A, we 
multiply the last inequality in (2.11) by A and we add it to the sum of the other two; 
the result is 

^2(1 -E) + cA{l + E)-{l + k){l + E) > , 
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hence + cA — {1 + k) > 0, which imphes A > X^. The second inequahty in (2.12) 
foUows, since by (2.11) cE > A{D-1 + E) > X+{D-1-E). This proves that V+ C 7^2- 

Finally, if 7(0;) = {a, P,w){x) is any trajectory of (2.1) in 7^2, then P'{x) < and 
a'{x) + DP'{x) < 0, hence P{x) and a{x) converge as a; ^ +00. Thus, proceeding as in 
the case D < 1, one shows that 7(0;) P2 as x ^ +00. This proves that 7^2 C 0,{P2). 
□ 

Lemma 2.2 ensures the existence of a nonnegative propagating front if the condi- 
tions (2.11) can be satisfied for some ^ > 0, £■ > 0. These conditions are equivalent 
to 

c > DA+l/A and c > A (^1 + 

for some E < {l + k)-'^{2DA^ - k). Therefore, given D>l,k>^, they can be fulfilled 
if and only if c > c{D, k), where 

c{D,k) = min max IdA+I/A , a(i+ t 11 • (2-13) 

In particular, setting A? = 1/D, we see that c{D, k) = 2VD if k < (3D - 1)/{SD - 2). 
Similarly, setting A'^ = k/D, we obtain c{D, k) < VD{Vk+l/Vk) for all D > 1, k > 1. 
Finally, straightforward calculations show that 

Inn c{D,k) = c*ik) , liin = c*(2/c) , 

where c*{k) is given by (1.8), and that lim c{D,k)/Vk exists for all D > 1. This 
proves the upper bounds in (1.9) for the critical speed c*{D,k). 



Remark. If /c < {SD - 1)/{SD - 2), the conditions (2.11) are fulfilled for all c > 2V-D 
if ^ = 1/y/D and E = (D — 1)/{2D — 1); this shows that the propagating front satisfies 
the bound 

(1-/?) <a< (1 - /5) (1 + ^^/^j • (2-14) 

To prove the lower bounds in (1.9), we use a similar argument. First, we verify as 
in the case D <1 that there exists no nonnegative propagating front if c < 2\fD. Thus, 
we assume that c > 2\fD^ and we define the region T^-s C Q+ by 

= {{a,P,w)\0< P <1 , 1-P<a , w< -B(5{1 - (5)] , (2.15) 
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for some B > 0. Then, if 

DB^ > 1 , DB^ -cB + 1 > , k- 2DB'^ > , (2.16) 

the vector field (2.1) on 97^3 is always directed into TZ^, except on the face f3 = 0. 
Indeed, this is easy to verify for the faces P = 1 and a = 1 — p. If w = —5/3(1 — /3), we 
have as in (2.8) 

w' + Bf3'{l - 2(3) = Bp{l -/3){^- B{1 - 2/3)) - ^a/3{l + fc/3) 

< -^/?(1 - p) {-{DB^ -cB + l) + p{2DB^ -k)) < , 

since a > 1 — (3. In addition, the conditions (2.16) ensure that the unstable manifold 
V+ is contained in TZ^ in a neighborhood of Pi, namely 

C + DX+ , ^ 

1 < r-^ , and B < X+ , 

c + A_|_ 

see (2.9). The first inequality is obvious since D > 1, and the second one follows by 
adding the last two inequalities in (2.16) : the result is DB^ + cB — {1 + k) < 0, hence 
B'^ + cB - {1 + k) < 0, which implies B < A+. 

According to these results, any trajectory on the unstable manifold V+ either re- 
mains in TZs for all a; e H or leaves TZ^ (and the positive sector Q+) by crossing the 
plane /? = 0. Now, the conditions (2.16) also imply that no trajectory of (2.1) can 
stay in TZs and converge to P2 as a; — > +00, except on the invariant line P = w = 
(which does not intersect V+.) Indeed, since the eigenvalues (2.4) satisfy A*o < A*± < 0, 
any trajectory in TZ^ \ {P = w = 0} converging to P2 becomes tangent to one of the 
eigenvectors w± as a; — > +00. In view of (2.5), (2.15), this is possible only if S < |//-|, 
in contradiction with the assumptions DB"^ > 1, DB"^ — cB + 1 > which imply 
B > IfJ,-]. Therefore, if (2.16) holds, the invariant manifold V+ necessarily crosses the 
plane P = and no nonnegative propagating front can exist. In particular, if A; > 2, we 
set B = ■sJk/{2D), and we conclude from (2.16) that no nonnegative propagating front 
exists if c < ^/D(^/kJ2 + y^2/k). This proves the lower bounds in (1.9) for the critical 
speed c*{D,k). 
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2.3. The case D > 1 : Existence of the Critical Speed 

Let D > 1. In this section, we show that the existence of a nonnegative propagating 
front for some value of the parameters c, k implies the same property for all c' > c, 
A;' < k. Thus we fix c > 0, A; > 0, and we assume that a, /? is a nonnegative solution of 
(1.2), (1.3). As in [BN], it is easy to verify that a'(x) > 0, P'{x) < and a(x)+/3{x) > 1 
for all a; e H. Setting w{x) = P'{x) as usual, we consider the bounded region TZ4 C Q+ 
delimited by the following four surfaces: 

Si = {{a, (3, w) \ w = 0} , S2 = {{ctj (3, w) \ a = 1 — (3} , 

53 = {{Xa{x),P{x),w{x))\x eJR, < A < 1} , 

54 = {(a(x),l3(x), iiw(x)) \ x eJR, 0</i<l}. 
We also note 7^4 = 7^4 \ {Pi}. 

Lemma 2.3. For all c' > c, < A;' < A;, the region TZ^ above is invariant under the Row 
of (2.1)', and V'+cTZaC r2'(P2). 

Remark. Here and in the sequel, (2.1)' denotes the vector field (2.1) with c,k replaced 
by c\k\ and similarly for V+ and n'(P2)- 

Proof. We proceed as in the proof of Lemma 2.1 or Lemma 2.2. First, it is straight- 
forward to verify that the vector field (2.1)' is directed into TZ4 on the surfaces and 
^2. If P e ^3 n a7^4, then P = {Xa{x), P{x),w{x)) for some x e H, A e [0, 1]. Using 
the equations (2.1) satisfied by a,P,w, it is easy to show that the vector 

N3{x,X) = (0 , D~^a{cw + al3{l + kl3)) , aw) , 

is normal to S3 at P and directed outside TZ^. On the other hand, for any c',A;' the 
vector field (2.1)' at P is given by 

Vi{x,X) = {c\l - Xa - P) - Dw , w , -D'^c^w + XaP{l + k'P))) . 

Taking the scalar product, we thus obtain 

N3{x, A) • Vd{x, A) = ^ ((c - c')w + (1 - X)ap + {k - Xk')aP^) < , 

since A < 1, c' > c and A;' < k, hence the vector field (2.1)' is directed into 7^4 on 
S3 n dTZ4. Similarly, if P = {a{x), P{x), iJ,w{x)) e ^4 for some a; e IR, /i e [0,1], an 
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exterior normal vector at P is given by N4^{x, ji) = (w^ , Dw'^ — cw(l — a — P) , 0), and 
the vector field (2.1)' at P reads 

Vi{x,ii) = {c\l - a- (3) - D/iw , tiw , -D-^{c'nw + a(3{l + k'(3))) . 

Therefore, we have N^i^x^fx) ■ Vi{x,n) = w'^{c' — cfx){l — a — (3) < 0, since /U < 1 and 
c' > c. This shows that TZ^ is invariant under the fiow of (2.1)' for all c' > < k. 

To prove that V+ C 7^4, we may clearly assume that either c' > c or /c' < k, since 
V+ C 7^.4 by construction. Then, as in (2.9), it is sufficient to verify that 

1 < J < + , and a; < A+ , 

c' + A+ c + A+ 

where A+ = A+(c', A;') is given by (2.2). The last inequality is satisfied because A+(c, k) 
is strictly decreasing in c and increasing in k, and the other relations follow since D > 1. 
Therefore, C 7^4 for all c' > c, k' < k. 

Finally, if 7(0:) = {a,^,w){x) is any trajectory of (2.1)' in 7^4, then ^'{x) < and 
a'{x) + DP'{x) < 0, hence /3{x) and a{x) converge as a; — > +00. Proceeding as in the 
previous cases, one shows that 7(0;) — > P2 as x — > +00. This proves that 7^4 C f2'(P2)- 
□ 

Using this result, we are now able to complete the proof of Theorem 1.1. 

Proof of Theorem 1.1. For any D > 1, A; > 0, let /(-D, k) be the set of values of 
c > 2\/D for which there exists a nonnegative solution of (1.2), (1.3). It is not difficult 
to verify that this set is closed, hence by Lemma 2.3 I{D,k) — [c*,-\-oo) for some 
c*{D,k) > 2Vd. It follows also from Lemma 2.3 that I{D,k) C I{D,k') if k > k\ 
hence the minimal speed c*{D,k) is a nondecreasing function of /c > 0. Finally, the 
upper and lower bounds (1.9) for c*{D,k) have been established in Section 2.2. This 
concludes the proof of Theorem 1.1. □ 
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3. Stability of the Propagating Fronts 



Throughout this section, we assume that D > 0, A; e [0, 1], c > 2-\/^, and we denote 



by a,P the solution of (1.2), (1.3) whose existence is ensured by Theorem 1.1 or Theo- 
rem 1.3. To prove Theorem 1.4, we shall control the behavior of the solutions of (1.12) 
in the function space Yg defined by the norm (1.13), (1.14). We begin with a standard 
local existence result: 

Lemma 3.1. Let D > 0, k e [0,1], c> 2-/D, and let (/o, go) e Y^. Then there exists a 
time ti>0 such that (1.12) has a unique solution (/, g) e C°([0, ti],Y^)nC^{{0, ti],Y^) 



is the generator of an analytic semigroup in Y^. Moreover, the nonlinearity : Y^ —> Ys 
in (1.12) is locally Lipschitz, uniformly on any bounded subset of Y^. Therefore, by 
Theorem 6.3.1 in [Pa], there exists a time ti > such that Eq.(1.12) has a unique 
classical solution (/, g) e C°([0, ti],Y^) n C^((0, ti],Y^) with initial data (/o, ^o)- □ 

Remark. In addition, the proof shows that, for any bounded subset B C Y^, the 
existence time ti > is bounded away from zero uniformly for all {fo,go) G B. It 
follows that the solution (/, g) either exists for all t e IR+ or leaves any bounded subset 
of Yg in finite time. 

In the sequel, we fix do > 0, ko > sufficiently small, and we assume that D G 
[1 — dojl + do], k e [0, ko], c> 2\/D. For e > sufficiently small, we make the following 
assumption: 

Hypothesis He- There exists a classical solution {f,g) of (1.12) defined on some time 
interval [0, T] and satisfying 



for all t e [0,T]. 

Under this assumption, we shall study the time evolution of some energy functionals 
which control the size of the solution (/, g) in Y^. In particular, if e, do, ko are sufficiently 
small, we shall show that the norm ||(/, 5f)||iA2 remains bounded on any time interval 



satisfying {f{0),g{0)) = {fo,go). 



Proof. It is straightforward to verify that the linear operator 




\\{m,gmY^ < 6, 
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[0, T] by a quantity depending only on the initial data. This result will be the main 
step in the proof of Theorem 1.4. As in [KR], [GR], it is convenient here to split the 
problem in two parts: First, we shall construct weighted functionals, with weight e^^, 
which control the perturbations (/, g) ahead of the propagating front. Then, we shall 
introduce unweighted functionals to describe the behavior of {f,g) behind the front. 



3.1. Weighted Functionals 

Let p{x) = e~*^, where s is given by (1-14). If (/, ^) is any solution of (1.12) in Y^, we 
define the weighted functions F, G, H by 

F{x,t) = p{x)-^f{x,t) , G{x,t) = p{x)-^g{x,t) , H = {l + ds^)F + G , (3.1) 

where d = D — 1. Then F,G,H e H^(IR), and a direct calculation shows that G,H 
satisfy the system 

dtG = (1 + d)dlG + iid^G - (7 + P)G + PH + pN{G, H) , 
dtH = dlH + (// + 2ds)d:,H - (1 + ds'^)H + d{dlG - 2sd^G) (3.2) 
+ ds^ (^(7 + P)G -PH- pN{G, H)^ , 



where 



II = - 4D , P = ^.[^ '^^11 , 7 = 1 - a(l + 2kp) . (3.3) 



(l + ds2) ' 
The nonlinearity N in (3.2) is defined by 



N{G, H) = p-^N ( P^^, PG ) . (3.4) 



Remarks. 

1. Both functions (3, 7 in (3.3) are close to (3 if \d\ and k are sufficiently small. Indeed, 
since < l/i:> = 1/(1 + d) by (1.14), we have 1 + ds^ = 1 + 0(1^1), hence 

PI (5 = l + 0{\d\+k) . (3.5) 

On the other hand, using the bounds (2.10), (2.14), it is straightforward to verify that 

(l-M|-2/c)/5<7</5(l + Y^) , (3.6) 



- 19 - 



if d > —1/2, hence 7//3 = 1 + 0{\d\ + k). In the sequel, we shaU always assume that 
\d\, k are sufficiently small so that /3 > 0, 7 > for all a; e M. 

2. In the limit x — > +00, the equations (3.2) reduce to 

dtG = {I + d)dlG + iid^G , 

dtH = dlH + (// + 2ds)da:H - (1 + ds^)H + d{dlG - 2sda,G) . 

For this limiting system, the "energy" J {G'^+H'^)dx is non-increasing in time if d = D—1 
is sufficiently small. Indeed, the diagonal term —(1 + ds^)H has a good sign, and the 
only cross term d{d'^G — 2sdxG) is a derivative and is multiplied by the small parameter 
d. This almost diagonal form in the limit x — > +00 explains our choice of the variables 
G, H instead of F, G. 

3. In the sequel, we shall use the estimate 

which holds for all c > 2\/D if \d\ and k are sufficiently small. To prove (3.7), we first 
note that I3'{x) + sl3{x) > for all x G IR. Indeed, using (1.2), (1.14), we see that the 
function z — (3' + s(3 satisfies the equation 

Dz' + (c - Ds)z + {a + kaf3 - 1)/? = . 

On the other hand, the bounds (2.10), (2.14) imply that a + ka^ - 1 < /3(-l + k+ \d\), 
hence Dz' + (c - Ds)z > for all x G IR if /c + \d\ < 1. Since z{-oo) = s > 0, 
this differential inequality implies that z{x) — (3'{x) + s(3{x) > for all x G IR, hence 
P{x)/p{x) is an increasing function. Now, recall that in the introduction we used the 
translation invariance of the problem to impose that P{x) — a{x) > 3/4 for all x < 0. 
In particular, we have /3{0) > 3/4, hence /3{x)/p{x) > /3{0) > 3/4 if a; > 0. Since is a 
decreasing function, we also have P{x) > /3(0) > 3/4 for all a; < 0. This proves (3.7). 

To control the evolution of G, H in H^(IR), we introduce the energy functionals 
Eo{t) = \ I {G^ + H^)dx, 

Ei(t) = ^ y (dG"2 + (7 + /3)G'2 + + (1 + ds')E''^ dx , 
where G' = d^G, H' = da^H. 
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Lemma 3.2. There exist do > 0, ko > and e > such that, if the hypothesis Tie 
above is satisfied, then there exists Ki > such that, for all t e (0, T], 

EQ{t) < -KiE,{t) . (3.8) 



Remarks. 

1. The constant Ki in (3.8) is independent of d e [—do, do], k e [0, /cq], c > 2y/D and 
T > 0. 

2. Here and in the sequel, we denote by " the time derivative to distinguish it from the 
space derivative ' . Unless stated otherwise, all the integrals are taken over the whole 
real line IR. 

Proof. Since {f,g) e C^{[0,T],Y,^)nC\{0,T],Y^), we have Eo e C^{[0,T])r\C\{0,T]). 
Using the first equation in (3.2) and integrating by parts, we obtain 



y = y (|-(1 + d)G'^ - (7 + P)G^ + PGH + pN{G, H)G^ dx 

< J (^-{1 + d)G'^ - (7 + /3/2)G2 + + pjv(G, H)G^ dx 

Similarly, using the second equation in (3.2) and integrating by parts, we have 

I H'^dx = ! - (1 + rfs2)iy2 - dG'H' - 2dsG'H) dx 

2 dt J 

+ ds^ j (yi + ^)GH - PH"^ - pN{G, H)H^ dx . 



(3.9) 



Since -dG'H' < \d\{G''^ + H'^)/2, -2dsG'H < \d\s{G'^ + H^), ds^GH < \d\s^{G^ + 
H'^)/2, and since < {1 + d)~^, there exists Co > such that 



1 d 

2Jtj^ 



J H^dx < - J {H'^ + H^) dx - ds^ j pN{G, H)H dx 

+ Co\d\J [h'^ + H^+ G"2 + (7 + p){G^ + H^)^ dx . 



(3.10) 



To bound the nonlinear terms in (3.9), (3.10), we first note that, due to the hy- 
pothesis 'Hei there exists Ci > such that 



p{x)\G{x,t)\ < CiePix) , 



(3.11) 
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for a\lxeM,te [0, T]. Indeed, iix <0, using (3.1), (3.7) and the embedding of Hi(IR) 
into L°°(1R), we find 

p{xf\G{x)\' = \g{x)\' < \\g\\l < WghWh < IbWh < IpixrMh • 

If a; > 0, we have by (3.7) 

1 fi 1 fi 1 fi 

p{xr\G{x)\' < -P{xr\G{x)\' < -p{xr\\G\\l < -PixfWGUG'h . 

Since \G'{x)\ < {s\g{x)\ + |^'(x)|)e"^, we have ||G||2||G"||2 < C||^|||.^ for some C > 0. 
Therefore, there exists Ci > such that p{x)\G{x)\ < Ci(3{x)\\g\\Y^ for aU x G M, and 
(3.11) foUows from the hypothesis Ti^. 

Now, using (3.4), (1.12), we have 

J p\N{G, H)G\ dx < J p\G\ + + kp\G\) + kaG^^ dx . 

Since \H - G\\G\ < {3G^ + H^)/2 and (1 + ds^)'^ = 1 + 0{\d\), it foUows from (3.11) 
that there exists C2 > such that 

/ p\N{G,H)G\dx < Cie [ {3 (1 + 2fe + kCie) + kG^ dx 

J J \2{1 + ds^) J ^3_^2) 

< Cse j p{G'^ + H^) dx . 

In a similar way, we obtain 

J p\N{G, H)H\ dx < Cae J (3{G^ + H'^) dx . (3.13) 

Therefore, combining (3.9), (3.10), (3.12), (3.13), and using (3.5), (3.6), we see that 
there exists C3 > such that 

Eo <- j (g''' + (7 + P/2)G^ + H'^ + ^H^^ dx 

+ Cs{\d\ + k + e) J (g"2 + (7 + p)G^ + H'^ + H^^ dx . 

In particular, if \d\, k and e are sufficiently small, there exists Ki > (independent of 
c > 2/D and T > 0) such that 

Eo < f (DG"' + (7 + /3)G' + iy" + (l + rfs2)iy') dx = -K^E^ . 
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This concludes the proof of Lemma 3.2. □ 

Lemma 3.3. There exist do > 0, ko > and e > such that, if the hypothesis He is 
satisfied, then there exists K2{c) > such that, for aU t G (0,T], 

Ei{t) < K2{c)Er{t) . (3.14) 

Remark. The constant i^2(c) is independent of T > 0, and behaves hke for large 
values of c, uniformly in (i G [—do, do] and k G [0, /cq]- 

Proof. We start from the identity 

h = \j^j ((1 + rf)G"' + (7 + /3)G'') cZa; = j (-(1 + d)^" + (7 + /3)G') Gda; . 

Using (3.2), we replace -(1 + (Z)G'" + (7 + ;5)G' with -G + fxG' + PH + pN{G, H) in the 
right-hand side. We obtain 

h = J (-G^ + i^G'G + (3HG + pN{G, H)G^ dx . 
Since iJiG'G < n^G'^ + G^/4 and ^HG < (3{H^ + G^)/2, we have 

h < J (^~G^ +p^G'^ + ^^{H^ +G^)+pN{G,H)G^ dx . (3.15) 

Similarly, we have the identity 

/2 = 11 y (iy'2 ^ ^ ds^)H^) dx = J {-H" + (1 + ds^)H)Hdx . 

Replacing —H" + (1 + ds^)H in the right-hand side with its expression obtained from 
the second equation in (3.2), we find 

h = J (--ff^ + (jw + '2ds)H'H + dG"H - 2dsG'H^ dx 

+ ds^ j (^(7 + /3)G - ^ff - pN{G, H)^ H dx . 

In view of (3.2), we also have 

G" = -^^^ (g - pG' + (7 + /§)G -PH- pN{G, i?)) . (3.16) 
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Replacing (3.16) into the expression of I2, we find 

h = J i^-H'^ + (/i + 2ds)H'H - 2dsG'H + ds^{^ + ^)GH - ds^^HH^ dx 

- (^Y^+ds^^ J pN{G,H)Hdx . 

Since {p + 2ds)H'H < H"^ /2 + {y.^ + 4d^s^)H'^ and \nG'H\ < n^G'^/2 + H^/2, it is 
easy to verify that there exists C4 > such that 



/2 < y (^-^H^ + n^H'^^ dx + C4\d\ J p\N{G,H)H\dx 
+ C4\d\ J [h'^ + + + + (1 + /i2)G"2 + (7 + 



(3.18) 



da; 



To bound the nonhnear terms in (3.15), (3.18), we proceed as in the proof of 
Lemma 3.2. We obtain 



j p\N{G, H)G\ dx < Cse J P{G^ + G^ + H^) dx , 
J p\N{G, H)H\ dx < Cse j p{G'^ + H'^ + H^) dx . 



(3.19) 



Therefore, combining (3.15), (3.18), (3.19) and using (3.5), (3.6), we see that there exists 
C5 > such that 



+ CM + k + e) J (i?" + H'^ + H^ + G^ + {1 + fx^)G'^ + (7 + p)G^) 



dx 



In particular, if \d\, k and e are sufficiently small, there exists Cq > such that Ei < 
Ce{l + fi^)Ei, hence (3.8) holds with -^"2(0) = Ce{l + /x^). This concludes the proof of 
Lemma 3.3. □ 
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3.2. Unweighted Functionals 

To control the perturbation (/, g) behind the front, we define h = f+g, and we consider 
the equations satisfied by /, h. From (1.12), we obtain 

dtf = dlf + cdj -df- a{l + 2kp)h - N{f, h - f) , 
dth = Ddlh + cd^h - ddlf , 

where 5 = (3 — a -\- k/3{/3 — 2a). As in the weighted case, the variables /, h have been 
chosen so that the system (3.20) becomes almost diagonal in the limit x — > — oo. To 
control the evolution of f,h in H^(IR), we define the functionals 

E2it) = ]- I {f + h^ + 2a^h'')dx, 
2 Jm 

Esit) = \ [ (r + O + f / {f + 2a'h')dx, 

where K > 6 is an arbitrary constant. 

Before computing the time derivative of E2,Es, we note that the additional term 
J a^h^dx in E2 satisfies 



d 
di 



j a^h^dx = -2D j a'^h'^dx + 2 J ah?dx-2d J a'^f'hdx , (3.21) 



where a{x) = {D{a'^)" — c{a^)')/2. In the sequel, we shall use the following two prop- 
erties of cr: 



i) Since a{x) ~ 6"^+^ as a; ^ —00, where A+ is given by (2.2), we have 

lim = 2DX%. - cA+ = {2D + -1-k 

' — a-'^x) 



a;— > — 00 



< {2D + 1){VD + l + k- \fDf -1-k 
= -2V2{V2 - if + 0{\d\ + k) , 

where we used the fact that A_|_ < y/D + 1 + k — ^/D for all c > 2\^. Therefore, if 
\d\ and k are sufficiently small, there exists G IR such that a{x) < — a^(x)/3 for all 
X < xq. Using the translation invariance of the problem, we may assume (without loss 
of generality) that xq > 0, hence 

(j{x) < ^ , for all x < . (3.22) 
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This condition is obviously compatible with the previous requirement that (3{x) — a{x) > 
3/4 when a; < 0. 

ii) For all a; e IR, we have 

a(x) < 2D{l + k)a^(x) . (3.23) 

Indeed, let z{x) = a'{x) — Vl + ka{x). Then z{x) ~ e'^+^(A+ — Vl + k) as a; — > — oo. 
Since A+ < Vl + k by (2.2), it follows that z{x) < if a; is sufficiently negative. On the 
other hand, in view of (1.2), the function z satisfies the differential inequality 



+ (c + Vl + k)z = a" + col - (1 + k)a - cVl + k a 



= -a{l - (3){l + k + k/3) - cVl + ka < , 

which implies that z{x) stays negative for all x eJR, hence < a'{x) < Vl + ka(x) for 
all a; e H. Using (1.2) again, we conclude 

a = Da''^ + Daa" - caa' = Da'^ + Da'^(3{l + k(3) - c(l + D)aa' < 2D{1 + k)a'^ , 
which proves (3.23). 

Now, we shall control the time evolution of E2 and £'3. 

Lemma 3.4. There exist do > 0, ko > and e > such that, if the hypothesis Tie is 
satisfied, then there exist > and K4 > such that, for all t e (0, T], 

Mt) < -KsE^it) + K^E^it) . (3.24) 

Remark. The constants and K4 are independent of d e [—do, do], k e [0,A;o], 
c > 2Vd and T > 0. 

Proof. Using (3.20) and integrating by parts, we obtain 
1 d 



2dt 



J fdx = I i-f" - 5f - a{l + 2kp)fh - N{f, h - f)f) dx 



< 



I - Sf + (^ + kp){f + a'h') - N{f, h - f)f^ dx . 



(3.25) 



Since — (5 = a — /3 + k(3{2a — (3) < a — P + ka^ < a — (3 + k, we have — 5 < 1 + /c for all 
a; e IR and, due to our choice of the origin, — 5 < — 3/4 + /c for all a; < 0. Therefore, 



l + k(3-5]fdx< [ (-\+2k^fdx+[ (I + 2k) fdx 



1 \ 7 ^°° 

- +2k] fdx+ - fdx 
IR \ 4 J 4 Jo 



(3.26) 
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To bound the nonlinear term in (3.25), we observe that 

mf,h-f)f\ = \f(h-f)il + 2kp) + kafih-ff + kfih-ff\ 

< \\h-f\\^{fil + 2k + k\\h-f\U + ka\f{h-f)\) . 

Since a\fih - /)| < a/2 + (p + a^h?)/2 and since \\h - f\\^ = \\g\\^ < \\g\\Y^ < e by 
the hypothesis He, there exists Ci > such that 



\N{f, h - f)f\ dx < Cie j if + a^h^) dx . (3.27) 
Combining (3.25), (3.26), (3.27), we obtain 
\jj f"dx< j [-r'-\p + \a^h''^dx + {C,e + 2k) J{f + a'h')dx 



1 



(3.28) 



+ -/ rdx. 







d 
dt 



Next, using (3.20) and integrating by parts, we find 
^±Jh'dx = J i-Dh" + df'h') dx< J (^{-D + M)/,'2 ^ dx . (3.29) 

Finally, integrating by parts in (3.21), we obtain 

J a^h^ dx = j {-2Da'^h'^ + 2(7h^ + 2da^f'h' + 2d{a^)'f'h) dx 

< j {{-2D + \d\)a'h'^ + {2a + \d\{a')')h'' + ^{0" + (a')')/") dx . 

(3.30) 

In view of (3.22), (3.23), we have 

2 / ah'^dx < a'^h'^ dx + 4:D{l + k) / a'^h'^ dx 

JjR 3 J-00 Jo ^2 g-j^-j 

~~^/ ^^^^dx+ (^^+4:D{l + k)^ a^h^dx. 

Replacing (3.31) into (3.30) and recalling that (ck^)' = 2aa' < 2\J\ + k a^, we thus find 

-f- / a^h^dx <- [ (2a^h''^ + -a'^hA dx + 5 [ a^h^ dx 

dtj J \ S J Jo (3 32) 

+ C2{\d\ + k) J a\h'^ + h^ + f'^) dx , 
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for some C2 > 0. 

Therefore, combining (3.28), (3.29), (3.32), we see that there exists C3 > such 



that 



E2<-J + \f + + (1 + 2a^)h'^^ dx 

+ C^{\d\ + k + e) J {f'^ + f + h''^ + a'^h^)dx + 5j^{f + a'^h^)dx 



(3.33) 



In particular, if \d\, k and e are sufficiently small, there exists K3 > (depending on 
K) such that 

E2 < + 5 / {f + h'')dx . 

Jo 

It remains to show that J^{P + h'^) dx < CEi for some C > 0. Using (3.1), (3.7), we 
have for all a; > 

f + h"^ < Hf + g"^) = SpHf'^ + G'^) < 3p(F2 + G'2) < 4/3(f2 + G'2) . 

Thus, if \d\ and k are sufficiently small, it follows from (3.1), (3.5) that there exist 
C4 > 0, C5 > such that P + h^ < C^PiG^ + H^) < 0^00^ + (1 + ds^)H^), for all 
X >0, hence J^{f^ + h?) dx < 2C5E1. This concludes the proof of Lemma 3.4. □ 

Lemma 3.5. There exist do > 0, ko > and e > such that, if the hypothesis Tie is 
satisfied, then for aiJ t e (0,T], 

Ez{t) < KK^Ei{t) . (3.34) 
Proof. Using (3.20) and integrating by parts, we obtain 

\jj f^dx = J + sfr + a{i + 2kp)hr + N{f, h - f)r) dx . 

Since |5| < 1 + k, we have dff" < f"^/4 + (1 + kfp and a{l + 2k(3)hf" < f"^/4 + 
a^{l + 2kYh? . Moreover, as in (3.27), we have 

j \Nif, h - f)f"\ dx < Ci6 j (/"^ + f + a'h') dx . 

Therefore, there exists Cq > such that 
1 d 



2dt 



J f'^ dx < J (^-\f"^ + f + a^h^^ dx + CQ{e + k) J (/"^ + f + a^h^) dx . 

(3.35) 



1^ 

2dt 
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Similarly, using (3.20) and integrating by parts, we find 

Combining these results, we see that, if k and e are sufficiently small, then 

^^y"(/" + /i")rfx < (l + C6(e + /c)) J{f + a%^)dx. (3.36) 

On the other hand, using (3.28) and (3.32), we obtain as in (3.33) 
^ ''{f + 2a^h^) dx < -K f ( /'2 + 1/2 + la2/,2 + 2a^h'A dx 



2 dt .n-" ' - . V 4" 6 



/pOO 
(/'^ + f + a'^h'^ + a^h^) dx + 5K J {f + a^h^) dx , 

for some Cy > 0. Since K > 6 by assumption, it follows that, if \d\, k and e are 
sufficiently small, then 

E^=\jJ if" + h") d^ + ^JtJ(f + «'^') < ^Kj\f + a-'h') dx . 

Proceeding as in the proof of Lemma 3.5, we thus obtain Ez < KK4E1. This concludes 
the proof of Lemma 3.5. □ 



3.3. Proof of Theorem 1.4. 

We first state two Corollaries which are direct consequences of the preceding Lemmas. 

Corollary 3.6. There exist do > 0, ko > and e > such that, if the hypothesis is 
satisfied, then there exists Ko{c) > independent ofT such that, for aU t G [0,T], 

\\{f{t),gm\Yi < i^o(c)||(/o,^o)|ln^ . (3.37) 

Proof. According to the four preceding Lemmas, we can choose do > 0, ko > and 
e > sufficiently small so that, if the hypothesis Ti^ is satisfied, then the differential 
inequalities (3.8), (3.14), (3.24) and (3.34) hold for t G [0,T]. In particular, the function 
E4 = Ki{Ei + E2 + E3) + [K2 + + KK4)Eq is non-increasing in time for t G [0, T]. 



- 29 - 



On the other hand, it is straightforward to verify that there exist K5 > and Kq > 
(independent of c, T) such that 

K4{f,g)\\l, <E,< K4{f,g)\\l, , 

hence (3.37) holds with Kq = {Ke/K^y/'^. Since K2 = 0{c^) as c ^ +00, the same is 
true for Kq, hence Kq = 0{c) as c — > +00. □ 

Corollary 3.7. There exist do > 0, ko > and e > such that, if the hypothesis He 
holds for all T > 0, then the solution {f,g) of (1-12) satisfies 

hm \\{d,f{t),d,gm\x^ = 0. (3.38) 

t— »+oo " 



Proof. We choose do > 0, ko > and e > as in Corollary 3.6. Using (3.8) and (3.14), 
we first note that the positive functions Eo{t) and KiEi{t) + K2Eo{t) are non-increasing 
in time, hence converge as t +00. In particular, Ei{t) has a nonnegative limit as 
t — > +00. Now, by (3.8), we have 



I 



+00 2 

Ei{t)dt = —{Eo{0) - Eo{+oo)) < +00 , 
-^1 

hence Ei{t) converges to as t — > +00. Similarly, combining (3.8), (3.24) and (3.34), 
we see that the positive functions KiE2{t) + K4Eo{t) and KiEs{t) + KK4^Eo{t) are 
non-increasing in time. Since £^o(^) converges, the same is true for E2(t) and E^{t). By 
(3.8) and (3.24), we have 



i 



-1-00 2 

E:,{t)dt < —— (K^iEoiO) - Eo{+oo)) + K^{E2{0) - E2{+oo))) < +00 , 
K1K3 



hence E3{t) converges to as t — > 00. Since \\{dxf,dxg)\W2 < K'j{Ei -\- E^) for some 
K7 > 0, (3.38) follows. ' □ 

We are now able to complete the proof of Theorem 1.4. 

Proof of Theorem 1.4. Let do > 0, /cq > 0, e > be as in the proof of Corollary 3.6, 
and let D e [1 - rfo,l + ^o], k e [0, fco], c > 2/D. We set eo = e/{2Ko), where 
Kq is given by Corollary 3.6. If (fo^go) e Ys satisfies || (/o, ^tq) ||Y-j2 < eo, then by 
Lemma 3.1 there exists a time ti > such that (1.12) has a unique classical solution 
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GC0([0,ti],i;2)nCH(0,ti],n2) satisfying (/(0),^(0)) = (/o,^o). By Corollary 3.6, 
it follows that 

\\{f{t),gm\Y^ < i^o||(/o,^o)||y,^ < 6/2 , (3.39) 

for all t G [0,ti]. Indeed, assume that there exists T G (0,ti] such that || (/(t), (/(t)) ||y2 < 
e for all t G [0,T) and || (/(T), (/(T)) ||y-2 = e. Then the hypothesis He is satisfied on 
[0, T], and Corollary 3.6 implies that (3.39) holds for t G [0, T], which is a contradiction. 
Therefore, || (/(t), (/(t)) ||y2 < e for all t G [0,ti], and (3.39) follows from Corollary 3.6. 
This shows that the solution {f,g) of (1.12) in satisfies (3.39) whenever it exists. 
Using the Remark after Lemma 3.1, we conclude that {f{t),g(t)) exists for all t > 
and satisfies (3.39). This proves (1.15), and (1-16) follows from Corollary 3.7. This 
concludes the proof of Theorem 1.4. □ 
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